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A simple teclinique for robust generation of vibrational Fock states in a chain of trapped ions 
is proposed. The method is fast and easy to implement, since only a single chirped laser pulse, 
simultaneously addressing all of the ions, is required. Furthermore, because the approach uses 
collective adiabatic passage, significant fluctuations in the intensity or frequency of the laser pulse 
can be tolerated, and the technique performs well even on the border of the Lamb-Dicke regime. We 
also demonstrate how this technique may be extended in order to create non-classical superposition 
states of the ions' collective motion and Greenberger-Horne-Zeilinger states of their internal states. 
Because only a single laser pulse is required, heating effects arising under realistic experimental 
conditions are negligibly small. 
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I. INTRODUCTION 

The interaction of a single quantum harmonic oscilla- 
tor with one or more two-level systems is a fundamen- 
tally important model in quantum physics ^ and offers 
an ideal arena for investigating the principles of superpo- 
sition and entanglement in mesoscopic systems 0. Many 
physical systems are well approximated by such a model, 
including, for example, atom-cavity systems chains 
of trapped ions with coupled internal and vibrational 
degrees of freedom quantum dots interacting with 
photonic crystal fibres or nanocavities Q and Joseph- 
son junctions interacting with microwave- resonators [6|. 
An area of active interest is to engineer non-trivial quan- 
tum states of the harmonic oscillator through a controlled 
manipulation of the two-level systems with which it in- 
teracts. This approach is well-established within cavity 
quantum electrodynamics where atoms passing through a 
cavity have been used to generate Schrodinger-cat states 
and Fock states [1| of the cavity field. 

The motional state of one or more trapped ions is an- 
other highly promising system in which to investigate 
non-classical harmonic oscillator states. When a chain of 
ions is cooled sufficiently, each vibrational normal mode 
closely approximates an ideal harmonic oscillator. By 
using appropriately tuned lasers, the transitions in the 
internal state of each ion are accompanied by changes in 
the vibrational state of a particular normal mode [J, Q . 
Furthermore, each mode couples very weakly to its en- 
vironment (via fluctuating charges on the trapping elec- 
trodes), so decoherence rates are typically much lower 
than in cavity QED systems. 

Many theoretical proposals exist for the creation of 
non-classical harmonic oscillator states in the motion 
of trapped-ions. In particular, the preparation of mo- 
tional Fock states (eigenstates of the harmonic oscil- 



lator) has been prop osed via the observation of quan- 
tum jumps pdl . Il2|. utilising adiabatic passage tech- 
niques [l3|, or using the method of trapping states 
[li, [3. Several non-classical states have been pre- 
pared experimentally and fully characterized by quan- 
tum state tomography, including Fock states, superposi- 
tions of Fock states. Schrodinger cat states and squeezed 
vacuum states [l^ [O, [H, Furthermore, the deco- 
herence of these states has been investigated by using 
specially engineered reservoirs [13, HH- However, despite 
several promising theoretical proposals for their genera- 
tion, to date, motional Fock states have only been suc- 
cessfully created by applying a sequence of Rabi 7r-pulses 
- tuned to the blue- and red-sidebands for a single ion 
[l6l. [2^. While this approach is conceptually simple, it 
is also technically demanding, since an iV-phonon Fock 
state requires N separate pulses of precise area and fre- 
quency, and each of these must be applied using a tightly- 
focussed beam addressing a single ion. Intensity and fre- 
quency fluctuations result in a reduced fidelity. 

In this paper, we propose a novel and very simple 
adiabatic passage technique for the generation of Fock 
states, requiring only one common laser pulse address- 
ing the whole ion chain. This technique can readily be 
adapted to the creation of Greenberger-Horne-Zeilinger 
(GHZ) states and superpositions of phonon Fock states. 
The approach has several significant advantages over ex- 
isting proposals, namely: (i) it is simple to implement, 
since only a single laser pulse is required; (ii) the state 
preparation is very fast, since it occurs in a single step; 
(iii) because there is only a very limited time for decoher- 
ence, the overall fidelity can be high; (iv) our technique 
is naturally robust against fluctuations in the laser Rabi 
frequency and detuning, since it utilises adiabatic pas- 
sage; (v) the technique holds up even when the first-order 
Lamb-Dicke approximation begins to break down. 
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The paper is arranged as follows. In Sec. |TT1 we intro- 
duce a model Hamiltonian and Hilbert space for the case 
of N trapped ions interacting with a common laser pulse. 
It is shown that by a judicious choice of basis (calculated 
using the generalized Morris-Shore transformation [1^), 
the dynamics can be confined to an (N + l)-dimensional 
subspace of the overall (2^-dimcnsional) Hilbert-space. 
Furthermore, by sweeping the laser detuning through res- 
onance with a motional sideband, the ground state can 
be connected adiabatically to the highest energy state 
which has exactly N vibrational quanta. This is shown 
explicitly for the first blue-sideband in Sec. IIII Al and 
the first red-sideband in Sec. IIIIBi using a 'complex- 
sech' laser pulse shape. In Sec. IIV| it is shown that our 
technique may also be adapted in order to engineer GHZ 
states and superpositions of phonon Fock states of the 
ion chain. In Sec. |V]we estimate vibrational heating ef- 
fects relevant to an experimental implementation of our 
technique. Finally, in Sec. IVI[ we conclude our findings. 



II. DEFINITION OF THE PROBLEM 

We consider N identical two-state ions confined in a 
linear trap and interacting with a common laser field. We 
assume that the ions are initially cooled to their ground 
state of motion [2^ and that the trapping frequencies are 
suitably chosen so that radial modes remain unexcited. 
The laser frequency ijJL{t) is tuned near to the centre- 
of-mass blue-sideband resonance, LOL{t) = ijJq + v — S{t), 
where ujo is the Bohr frequency of the transition, i> is the 
trap frequency and S (t) is a time-dependent detuning. 
The chain of ions possesses N axial vibrational degrees of 
freedom and the corresponding vibrational normal mode 
frequencies are denoted Vp (where p = 1 . . . N and vi = 
v). If the laser field has equal intensity at the position 
of each ion - for example if it is directed along the trap 
axis - then, after applying the electric-dipole and optical 
rotating- wave approximations, the interaction between 
the ions and the laser field is described by the following 
Hamiltonian 0, [l^l : 



H/(t) 



icj). 



=^^|a+exp [^j\r)dr-^ut 

X exp 1^51 [ape^'""'' + ^e'""*]^ + /i.c.j. 

(1) 

Here a'^ = (Oj| and aj = |0j) (lj| are the rais- 
ing and lowering operators for the internal states of 
the j^^ ion, and Op are respectively the creation 
and annihilation operators for phonons in the p*^ mode. 
?7 = \J hk"^ cos^ d/2Mv is the single-ion Lamb-Dicke pa- 
rameter, with k being the laser wavenumber, 9 the angle 
between the trap axis and the direction of the laser beam 
and M the mass of the ion. is a dimensionless cou- 
pling constant that depends both on the particular ion 



and the mode under consideration. The function ri(t) is 
the (real-valued) time-dependent Rabi frequency and 

0j = 0^ - A:a;° cos 61 (2) 

is the phase of the laser pulse at the equilibrium position 
x° of the j*'' ion. 

The laser parameters are chosen to satisfy (5(t)/i/ <C 1, 
ri(i) < V and in this limit, off-resonant transitions to 
higher-order vibrational sidebands have a negligible effect 
on the dynamics. Therefore, to good approximation, we 
may write: 



a'^ exp 



X exp 



N 



\ae-""^ + a^e'"^] 



h.c. 



(3) 



where for notational convenience, 



a 



a 



; write a = ai, 

and we have used the fact that k-{ — 1/ \/N V j 
[9|. Equation Q simplifies considerably when 
fyV n + 1 <C 1; where n is the number of phonons in the 
ccntre-of-mass mode - a condition known as the Lamb- 
Dicke limit. After expanding in powers of 77 and applying 
the vibrational rotating-wave approximation, the domi- 
nant contribution to Eq. ^ is [2^ 



N 



h.c. 



,(4) 



where we have introduced a scaled coupling, 
g (t) = r]fl (t) between the internal and mo- 
tional degrees of freedom. We remark that although 
the following analysis uses the simplified Hamiltonian 
(|^. the numerical simulations presented below are 
based on the Hamiltonian given in Eq. ([3]). The effect 
of the additional counter-rotating terms and terms 
higher-order in rj is investigated in Sec. IIII CI and found 
to be negligible for a wide range of parameters. The 
above Hamiltonian describes the anti-Jaynes-Cummings 
model, which conserves the difference between the 
number of ionic and vibrational excitations - a feature 
characteristic of blue-sideband excitation. We assume 
that the initial state is IO1O2 . . . Oat) |0), where \n) is 
the Fock state with n phonons (n = 0,1,2...) and 
|0i02...0Ar) is the collective ground state of the ions. 
Each manifold with n phonons and n ionic excitations is 
C,^-fold degenerate, where = N\/[n\{N - n)!]. An 
example linkage pattern for three ions is illustrated in 

Fig.m 

The entire Hilbert space Ti. is spanned by the following 
basis states: 



|0i02...0w)|0), 

|Oi...lfe...OAr)|l), 
|Ol...lfe...l„...OAr) |0), 



(5a) 
(5b) 
(5c) 

(5d) 
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|111>|3> 



|111)|3) 



S5{t) 



|011)|2) 




|110>|2> / |101>|2> 



|100>|1> 



1000) |0> 

FIG. 1: (Color online.) Linkage pattern of the collective states 
of three trapped ions coupled to their common centre-of-mass 
mode. The difference between the number of ionic excitations 
and the number of phonons is conserved. The laser beam is 
tuned near to the blue-sideband resonance with detuning S (t). 
The coupling between the manifolds indexed by n and n + 1 
is g (t) ^/nTfT. 

and has dimension diniTi = J2n=o^n = 2^. We now 
perform a time-dependent phase transformation 



U(t) = cxp 



N 



5(r)dr -0,+7r/2 1 1,) (1,1 



(6) 

and re-express the Hamiltonian Q as 

H, (t) = U(t)H, (t) Ut(i) + ih^^V^t), (7) 

so that in the basis defined in ^ , H/ (t) has a tridiagonal 
block-matrix form 



H/ {t) = h 



Do 


Vo.i 


... 








vS,t 


Di 


Vi,2 ... 











1,2 


D2 ... 









D 



V 



N-l.N 

(8) 

dimensional intcrac- 
and (n + l)-phonon 



Here V„.„+i (t) is a (C^ x C„+i 
tion matrix which couples the n 
manifolds. (N—n) of the elements in each row and (n+l) 
of the elements in each column of V„^„_|_i (t) are non-zero, 
and from the symmetry of Eq. (|4]), it can be seen that all 
of these non-zero elements are equal to g(t)\/n + 1. D„ 
is a C^-dimensional diagonal matrix whose elements are 
equal to the respective shared detuning of the n-phonon 
manifold 



3S{t) 



25{t) 




^2,3(t) 



Al,2(t) 



Ao,i(t) 



|000>|0> 

FIG. 2: (Color online.) Example of the Morris-Shore trans- 
formation for three ions. For the generation of Fock states, 
the dynamics is confined to a four-state chain, where the de- 
tuning of level n is n5{t) and the coupling between levels n 
and (n + l) is A„,„+i(f). 



where 1„ is the unit matrix. 

On first inspection, an analytic treatment of the dy- 
namics may appear to be highly demanding for more 
than a few ions, since in the original basis (O, the dy- 
namics is described by 2^ coupled differential equations. 
Fortunately however, by utilising the multilevel Morris- 
Shore transformation |23|, it is possible to decouple this 
vast Hilbert space into a set of smaller closed subspaces, 
which simplifies the problem dramatically. Such a mul- 
tilevel Morris-Shore factorisation exists if and only if the 
following commutation relations are satisfied for all pairs 
of neighbouring coupling-blocks [n = 1, . . . , N — 1) 



V 



n—l,n 



n V n,n+l v n,n+l 



0. 



(10) 



It was recently shown that these commutation relations 
are automatically satisfied by the Hamiltonian given in 
Eq. ^ [26|; the 2^-dimensional Hilbert space can there- 
fore be factorized into a set of independent chainwise link- 
ages, as shown in Fig. O Furthermore, for an appropri- 
ately chosen initial state, the dynamics is wholly confined 
within the largest of these subspaces - an {N + l)-level 
ladder. The states which make up this chain are equal to 



U(0|Vt^,f)|n) (n: 



(11) 



where the symmetric Dicke states 
follows f27l,|28l 




0,...,iV), 
\W^) are defined as 



(12) 



D„ (t) = nS (t) 1, 



(9) 



and {Pfe} above denotes all possible distinct permuta- 
tions of the n ionic excitations. A full derivation of the 
Morris-Shore factorisation described above can be found 
in [2^, but as a further justification, we note that H.i{t) 
is symmetric with respect to the interchange of any two 
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ions. Therefore, states which are unchanged by a permu- 
tation of the ions can only be coupled by the Hamiltonian 
([HI) to other states with this same symmetry. The initial 
and final states are both symmetric product states and 
therefore, the system begins and remains in the subspace 
defined by Eq. (fTTj) throughout the dynamics. 

When expressed in the factorized Morris-Shore ba- 
sis, the entire evolution of the system is thus confined 
to an {N + l)-level ladder, where the lowest state is 
|0i02...0Ar) |0), the highest is |lil2...1jv) |^), and aU 
intermediate levels (n = l,...,iV— 1) are the multipar- 
tite entangled states ([TT]) . Each of the basis states con- 
tributing to level n is connected to {N — n) of the basis 
states which make up level (n + l) in the ladder, and this 
coupling has strength g{t)^/n + 1. After a short calcula- 
tion, the coupling between adjacent levels in the Morris- 
Shore chain is found to be 



A„,„+i {t) ^ g {t) {n + 1) VN-n. 



(13) 



Within the closed subspace formed by the {N + l)-level 
Morris-Shore ladder, the Hamiltonian reduces to 



H 



N+l 



(<) = h 



Ao,i 
Ao.i S Ai.2 
Ai 2 25 



... {N -1)6 Xn-i.n 
... \n-in NS 

(14) 

As an alternative to the Morris-Shore formalism, we 
note that the state-space factorisation presented above 
can also be conveniently derived via a careful in- 
spection of the Hamiltonian. A convenient route 
is to rewrite ^ in terms of the combined ionic 
pseudo-spin operators, Jz = ^ SjLi ~ 

J+ = Ef=i Il.)(0,l, J~ = Ef=i I0,)(l,|, as in Refs. H, 
[sol . [3]| . Our approach complements the collective-spin 
picture, and is also more general, in that it is not strictly 
limited to the case of equal couplings [1^. 



III. CREATION OF MOTIONAL FOCK STATES 



A. Creation of Fock states via a level crossing on 
the blue-sideband 



1. Bow-tie ladder model 

We wish to transfer the population of the initial state 
|0i02...0Ar) |0) to state |lil2...lAr) \N). This can be 
done very efficiently by creating a 'bow-tie' level crossing 
linkage pattern, wherein all energies [the diagonal ele- 
ments of the Hamiltonian p4)) ] cross at the same instant 
of time, and enforcing adiabatic evolution. The conve- 
nience of the bow-tie linkage pattern for our multistate 
Morris-Shore ladder is determined by the existence of 



an adiabatic state - a time-dependent eigenstate of the 
Hamiltonian (jl4p - which connects the desired initial and 
final states I Oi02...0jv) |0) and |lil2 • . . Ijv) \N). The en- 
ergy of this adiabatic states is either below (for up-chirp, 
when the slope Sq > 0) or above (for down-chirp, when 
the slope Sq < 0) all other adiabatic energies. By im- 
posing adiabatic evolution conditions, we can force the 
system to follow this adiabatic state, thereby arriving in 
the desired Fock state |lil2 . . . l^v) |A^) in the end. 



Detuning 



adiabatic energy 



111- 




>-Time 
^ coupling, 



FIG. 3: (Color online.) The motional ground state and A^- 
phonon Fock state of the ion chain can be connected adia- 
batically using a level-crossing, providing that the couplings 
are chosen to satisfy the adiabaticity condition, (|17|l . To em- 
phasise the 'bow-tie' structure of the level-crossing, the en- 
ergies of all of the diabatic states (blue curves) are shifted 
by -N5{t)/2 from those in Eq. The red dotted curve 

shows the particular adiabatic pathway used for Fock state 
preparation. 

To be specific, we choose to use a hyperbolic-secant 
pulse shape and a hyperbolic-tangent detuning. 



n{t) = Qa scch{t/T), 
5{t) = Sotanh {t/T) . 



(15a) 
(15b) 



When only two states are involved, this model is exactly 
soluble analytically and is known as the complex-sech 
model in nuclear magnetic resonance |32| and the AUen- 
Eberly-Hioe model in quantum optics [33| ; these are spe- 
cial cases of the more general Demkov-Kunike model [38| . 
The model ([T5|) can therefore be viewed as an (N + 1)- 
state generalization of this model. 

There are several analytic solutions available in the 
literature for multistate bow-tie level-crossing problems. 
Analytic bow-tie solutions have been found for three [34 1 
and N states [H, [s^ 113] , in which the coupling is con- 
stant and the detunings are linear in time. However, 
in these solutions a single state is assumed to interact 
with all other states, which in turn do not interact with 
each other; this linkage is unfortunately different from 
the chainwise linkage in our Morris-Shore ladder. A 
chainwise-coupled bow-tie linkage pattern occurs in an 
analytic model of an rf-pulse controlled Bose-Einstein 
condensate output coupler [i^, |4l| ; here the underlying 
SU(2) symmetry allows for a reduction to an equivalent 
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soluble two-state problem. This symmetry, however, im- 
poses strict relations between the couplings, which un- 
fortunately are not met in our Morris-Shore ladder. 

Compared to other types of level crossings, for in- 
stance the Landau-Zener model [s^, the sech-tanh model 
(|15p produces a relatively fast adiabatic passage, with a 
very convenient and controllable excitation-profile, with 
a very high probability in a certain range around res- 
onance and rapidly vanishing probability outside this 
range; the high-probability range is efficiently controlled 
by the chirp-rate i5o [4^ . A flat and narrow excitation 
profile is particularly relevant here, since it is impor- 
tant that our laser pulse does not excite unwanted vi- 
brational sidebands, and this point is discussed further 
in Sec. [mCl 



2. Adiabatic conditions 

Because a multistate analytic solution for the model 
psp does not appear possible, we resort to numerical 
simulations; however, we shall make use of analytic in- 
sights derived from the above models in order to esti- 
mate the adiabaticity conditions. A transfer efficiency of 
1 — e (with e <C 1) in the exactly soluble two-state case 
requires the coupling strength and the detuning to satisfy 



21n(l/e) 



> 



ttSoT > ln(l/e). 



(16) 



where the signs of 6o and f^o are unimportant and as- 
sumed positive for definiteness. We have verified numer- 
ically that in the multistate Morris-Shore bow-tie ladder 
this condition is amended to 



2N ln(l/e) 



> 



7i-<5or>iVln(l/e). 



(17) 



The important conclusion from here is that both the 
chirp rate So and the peak Rabi frequency flo have to 
increase linearly with the number of states N for the fi- 
delity to remain above the value 1 — e. We note that the 
linear dependence of Qq on N agrees with that found in 
other systems involving multiple level crossings [43| and 
other types of multistate adiabatic passage [44| . 



3. Mechanism of Fock state creation 

The technique operates as follows. First, we apply 
a common laser pulse to all N ions, tuned near to the 
blue sideband. By adiabatically sweeping the detuning 
S (t) through resonance, all population will be transferred 
from the lowest to the highest energy state, [sij 



|0i02...0Ar)|0) 



I1112 



1 



\N) 



(18) 



In an optional second step, to create number-states 
of higher phonon number, the laser is tuned near the 
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FIG. 4: (Color online.) Time evolution of the populations 
Pq and Pg which correspond to the chain containing zero 
and eight motional quanta, for the case of = 8 ions. The 
Schrodinger equation was solved using the Hamiltonian given 
in Eq. i.e. neither the Lamb-Dicke approximation, nor 
the vibrational rotating- wave approximation were applied. 
The time-dependences of the Rabi frequency Q{t) and the 
detuning S{t) were of sech-tanh form (I15|l . and the parame- 
ters were: g{0)T = 12.5, 60T = 10, 77 = 0.3, g{0) = u/20tt. 
The final fidelity is 99.8%. 



carrier or red-sideband resonance, in order to return the 
population adiabatically, as follows: 

(a) for the carrier transition 

|lil2...1jv)|iV) |0i02...0jv)|iV); (19) 

(b) for the red-sideband transition 

|lll2...lAr)|7V)^|0i02...0Ar)|2Ar). (20) 

Hence, by repeatedly applying the blue-carrier or blue- 
red transition cycles, it is possible to create the phonon 
Fock states \kN) (with fc = 1,2, . . .). 

Figure 3] shows the time evolution of the populations 
Po and Ps of the states |0i . . . Os) |0) and |li . . . Is) |8) for 
the case of eight ions. The initial ground-state popula- 
tion is completely transferred through several intermedi- 
ate entangled states to the highest-energy state which is 
an eight-phonon Fock state. We note that r] = 0.3 for 
the example given in Fig. 4. Despite this relatively high 
value of the Lamb-Dicke parameter, the technique still 
works with a high fidelity. 

Since the technique described above uses adiabatic pas- 
sage, the transfer efficiency remains high even for consid- 
erable variations in the laser intensity and detuning - 
providing that the condition (fT7|) is satisfied. To demon- 
strate this fact. Fig. [5]shows a contour plot of the transfer 
efficiency as a function of SqT and g{Q)T. The plot has a 
broad island of high transfer efficiency in the top-right- 
hand corner, and in this region the process is insensitive 
to fluctuations in the experimental parameters. This ro- 
bustness represents an advantage over conventional tech- 
niques requiring precise pulse areas and frequencies. 
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FIG. 5: (Color online.) Numerically calculated transfer ef- 
ficiency for the creation of a five-phonon Fock state. The 
contour plot shows — log(l — J-), (where T is the fidelity of 
Fock-state creation), plotted versus the scaled coupling g{0)T 
and scaled detuning 5oT. The values of the Lamb-Dicke pa- 
rameter are respectively: (a) rj = 0.1; (b) rj = 0.5. The simu- 
lations solve the Schrodinger equation with Hamiltonian (|3} 
and neither the Lamb-Dicke approximation, nor the vibra- 
tional rotating-wave approximation is employed. 



B. Creation of Fock states via a level crossing on 
the red-sideband 



We note that the same Fock states can also be created 
in a single step using the red rather than blue motional 
sideband. In this case, the Hamiltonian is of Jaynes- 
Cummings, rather than anti-Jayncs-Cummings type 



H 



N 



-)dT — i(f>j -\-i7T /2 



h.c. 



(21) 



which conserves the total number of excitations. The 
gen eralized Morris-Shore transformation again applies 
[2q but in this case, the Morris-Shore states in Eq. pT|) 
are replaced by 



Uit)\W^^^)\n) (n = 0,...,iV). 



(22) 



Therefore an iV-phonon Fock state can be created as fol- 
lows. Firstly, the ions are initialized in their upper state, 
and the chain cooled to its ground state of motion, i.e. 
II1I2 . . . lAr)|0). Subsequently, a common laser pulse is 
applied, and swept adiabatically through resonance with 



the red sideband. This transfers the system sequentially 
through the ladder of intermediate entangled states given 
in Eq. ^ and into the target state, IO1O2 . . .Ojv>|^)- 
Higher-phonon states \kN) may be created by using the 
red-carrier or red-blue transition cycles in the manner 
outlined in subsection IIII Al 



C. Performance beyond the Lamb-Dicke limit 

It is clear from Figs. 4 and 5 that while our technique 
is designed to operate in the Lamb-Dicke regime, it can 
tolerate significant values rj and still perform with high 
fidelity. This represents another clear advantage of the 
adiabatic technique presented in this article over conven- 
tional approaches to motional state engineering. To see 
why the technique can function with a range of rj values, 
we note that the terms present in Eq. ^ which were dis- 
carded en-route to Eq. ^ can be categorised into three 
groups, as follows: (i) resonant terms involving only op- 
erators from the centre-of-mass mode; (ii) non-resonant 
terms involving only operators from the centre-of-mass 
mode; (iii) non-resonant terms involving operators corre- 
sponding to higher-order axial normal modes. 

The three groups of terms have the following effects. 
Terms in category (i) simply adjust the strength of the 
Morris-Shore couplings between different levels in the 
longest chain. It is well-known that adiabatic techniques 
can tolerate significant variations in coupling strength 
and so, unsurprisingly, these terms do not harm the pro- 
cedure. Category (ii) terms have a negligible effect pro- 
viding that the vibrational rotating-wave approximation 
holds. In fact, application of the rotating-wave approxi- 
mation is favoured both by a small Lamb-Dicke parame- 
ter, and also by the condition of adiabatic dynamics and 
these terms have a negligible effect even for flo ^ v/y/N 
in this regime. 

The terms in category (iii) correspond to multi-phonon 



transitions, at sideband frequencies (±fp ± + [52l |. 
where p and q run over all vibrational modes other than 
the centre-of-mass mode importantly, these are all 
non-resonant transitions. Because the scch-tanh pulse 
shape chosen in this paper is extremely selective, with 
an excitation profile that is both flat and narrow, ex- 
citation of off-resonant modes is strongly inhibited [i^l . 
In order to qualify this statement, we note that for the 
high-fidelity region in Fig. [SJi [i.e. g{0)T ~ 6, SqT - 10, 
g{Q) ~ i'/20tt], the maximum pulse detuning is on the 
level of 2.6% of the trap-frequency. By contrast, for a 
chain of ten ions, and taking into account the terms up 
to order 77^ , we find that the motional-sideband with fre- 
quency closest to the sideband we address is detuned by 
approximately 11% of the trap- frequency - far enough 
from resonance to remain unexcited, due to the high se- 
lectivity of the chosen sech-tanh pulse shape. Taking 
into account those terms proportional to 77^, the nearest 
motional sideband is still off-resonance by 4.7% of the 
trap frequency, and also remains unexcited. Each sue- 
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cessive group of terms also has a coupling strength that 
is of a factor of t] weaker again that the group before - 
further aiding selectivity. Therefore, from our numerical 
simulations, and also from the arguments outlined above, 
we conclude that our proposed technique should operate 
with high fidelity, even on the border of the Lamb-Dicke 
regime. 

IV. CREATION OF MOTIONAL FOCK STATE 
SUPERPOSITIONS AND INTERNAL GHZ 
STATES 

A. Creation of ionic GHZ states 

A modified version of the above tec hniq ue can be used 
to create generalized GHZ states [1^, |46| |. This is done 
in a two-stage process. First the number-state superpo- 
sition 

IO1O2 . . . 02„+i) (co |0) + ci \2n + 1)) (23) 

is prepared, and then the motional excitations are 
mapped onto the internal state of the ion chain. The 
state (f23|) can be created as follows. Initially we prepare 
the system in the product state IO1O2 . . . 02n+i) \n) us- 
ing the technique described above, and then perform a 
controlled rotation of the state of the first ion using the 
carrier-frequency transition, to obtain 

(co|0i)+ci|li))|02...02„+i)|n) (24) 

for even n or 

(ci|Oi)+co|li))|02...02„+i)|n) (25) 

for odd n. Next, a sequence of avoided level crossings 
is induced by addressing the first ion, with pulses alter- 
nating between the first red- and blue-sidebands (ending 
with the red sideband). We note that such level crossings 
have been performed experimentally by Wundcrlich et al, 
each with a fidelity above 99% [13] . During each avoided 
crossing, a phonon is removed or added, depending on 
the internal state of the first ion, and after {n + 1) cross- 
ings, the state ([^ is reached. This sequence of avoided 
level crossings is shown schematically in Fig. [51 

In the return step, all {2n + 1) ions are addressed si- 
multaneously by a laser beam swept through the red- 
sideband transition in the manner described in Sec. IIIIBI 
This maps the phonons onto ionic excitations and gives 

IO1O2 . . . 02„+i) (co |0) + ci \2n + 1)) ^ |GHZ) |0) , (26) 

where 

|GHZ) = Co IO1O2 . . . 02„+i) + cie'« II1I2 . . . l2„+i) (27) 

is a generalized GHZ state with arbitrary amplitudes, 
and the known phase ^ is specified in Eq. pO|) below. 
The above technique creates GHZ states of odd numbers 
of ions, from which even-numbered GHZ states can be 
derived using only single- ion rotations and measurement. 



fled Blue Red fled Blue 

sideband sideband sideband sideband sideband 




FIG. 6: (Color online.) Sequence of level crossings used to 
prepare the state (|23p for even n. The routes taken by the 
adiabatic states are shown as red dashed curves, while the 
bare energies are shown as solid blue curves. The initial state 
of the first ion is: (a) internal state |1), (b) internal state |0). 
After the n crossings shown above, a final level-crossing on 
the red sideband has the effect |l)|27i) ^ |0)|2n + l), while 
leaving the state |0)|0) unchanged. 



B. Creation of superpositions of motional Fock 
states 

From the preceding analysis, it is clear that our tech- 
nique can also be used to convert entanglement between 
the ions' internal states into a coherent superposition of 
two different states of their collective motion. Beginning 
in a generalized A'^-ion GHZ state and applying a common 
pulse to the whole chain (swept through the red-sideband 
transition) gives 

|GHZ) |0) ^ IO1O2 ... 0^) (co |0) + cie'€ |7V)) , (28) 

and applying the blue-carrier transition cycle k further 
times creates the following superposition state 

l^-fe) = IO1O2 ... 0^) (co \kN) + cie^«' |(fc + I) AT)) . 

(29) 



C. Phase factors 

In general, the phase factor appearing in Eqs. (|27p . 
(PS)) and has two components, 

e = ^ r E{t')dt' -J2^Ah\{\i'f)~M^ (30) 
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(b) 




"^"'^«;j||ll^)llltt0 



FIG. 7: (Color online.) Alternative experimental routes to the 
creation of vibrational number states using adiabatic passage: 
(a) Common addressing of all ions; (b) Individual addressing 
of each ion; (c) Repeated addressing of a specific ion in the 
chain; (d) Repeated addressing of a single trapped ion. Ap- 
proaches (a) and (d) require that the centre of mass mode is 
used, whereas (b) and (c) may be implemented using higher- 
order vibrational modes. 



where E{t) is the eigenvalue of the Hamiltonian cor- 
responding to the state used for the adiabatic transfer 
process, and (j>j is given in Eq. ([2]). In situations where 
the whole state accumulates a global phase, ^ has been 
omitted above, but in cases where only a constituent part 
of the total state acquires this phase (i.e. in Sec. lIVp . it 
has been included explicitly. 

The dynamical part of ^ [the integral in Eq. ([50)1 ] 
can be controlled by choosing appropriate values for the 
area and detuning of the laser pulse. Alternatively, if a 
Raman-coupled hyperfine qubit encoding is chosen, then 
the adiabatic passage described in section [IT] can be per- 
formed using a multi-ion dark state. In this case, E{t) is 
identically zero throughout the transfer process and the 
state acquires no dynamical phase. Full details of this 
dark-state technique will be described elsewhere [sO] . 

The second term in Eq. ([50)1 arises because the ions 
are not equally spaced along the trap axis. This means 
that using a single laser pulse, it is impossible to set all 
of the individual phase factors e*'^^ to unity. However, by 
symmetry, it is possible to set "/"i = 0- Hence when 

the whole chain is addressed simultaneously, this second 
term in (|30p always vanishes. For individual addressing 
of each ion, this term can always be chosen to be zero, 
even when only part of the chain is used. 



V. HEATING EFFECTS 

Vibrational heating is one of the major limiting fac- 
tors in ion-trap quantum information. Here we esti- 
mate the heating effects in the proposed global addressing 
technique and make a comparison with three alternative 
adiabatic-passage techniques: one that uses N ions and 
local addressing (which is a variation of our technique 
and allows us to use higher-order phonon modes), and 
two variations on the traditional one that uses a sin- 
gle ion. These different implementations arc sketched 
schematically in Fig. [71 



1. Center- of -mass mode and global addressing 

In writing Eq. (j4]), we have assumed that the coupling 
git) between the internal and motional states is the same 
for each ion, and produced by a single laser pulse, tuned 
close to the blue sideband of the centre-of-mass mode. 
While using a single laser pulse makes the technique ex- 
tremely simple to implement, the heating rates for the 
centre-of-mass mode are typically larger than for higher- 
order vibrational modes [2J, . 

As follows from the adiabatic condition (fTT)) . the re- 
quired pulse area A — nfloT of a sech-tanh pulse for 
fidelity 1 — e scales in proportion to N and ln(l/e). For 
a fixed peak Rabi frequency Qq, the characteristic pulse 
width for our global-addressing technique is therefore 



T = 



V2Nln{l/e) _ VN ln(l/e) 



TTr/flQ 



V2 ngiO) 



(31) 



The total interaction time for a sech pulse can be taken 
as Tgiobai « lOr. 

In order to estimate deleterious effects of heating on 
our scheme, we assume 17(0) ~ v jlQix. This is the rea- 
sonable maximal value of the coupling (implying minimal 
pulse duration), for which no extraneous excitations in 
other vibrational modes are created Q and the rotating- 
wave approximation is satisfied [sst . For a trap frequency 
of vjl-K = 4MHz as in [i^, we find that total interaction 
time of Tgiobai ~ 150 ^s is sufficient for 0.9999 fidelity 
(e ~ 10^"*) for N = 8 ions. With a heating rate per ion 
of approximately 5 phonon/sec [49| . we predict heating 
of about 6 X 10~^ phonons in total. 



2. Higher vibrational modes and local addressing 

Decoherence effects may be further reduced by using a 
higher-order vibrational mode with a lower heating rate, 
if each ion is addressed individually, i.e., 

N 



Each individual coupling strength, 

9 jit) = 



2VN 



(32) 



(33) 



now contains the dimcnsionlcss factor defined in Eq. 
(H]) [13] ■ By tuning the relative amplitude of each indi- 
vidual laser appropriately, it is possible to produce the 
Hamiltonian 



3. Using a single ion 

Alternatively, as a trade-off between the technical sim- 
plicity of using a single laser pulse and the low heating 



9 



rates afforded by higher-order modes, one can use a sin- 
gle ion repeatedly instead of N ions in a single step; then 
phonons can be added or subtracted sequentially instead 
of all in one step. In this case, the difference between 
addressing one out of a chain of N ions (Fig. [T);) and a 
single trapped ion (Fig. [7ji) is the heating rate, which 
scales linearly with N. Below, for clarity, we consider the 
approach involving one ion from a chain of N ions (Fig. 
^) , and compare this to the global addressing technique 
discussed in subsection lV II (Fig. [7^). 

The advantages of using a single ion are the lower heat- 
ing rates of higher-order modes, and the use of a single 
laser. It also permits the creation of number states of ar- 
bitrary phonon number, rather than in multiples of the 
total number of ions. However, we stress that the single- 
ion multiple-step scenario requires a significantly higher 
fidelity in each step, 1 — e/N instead of 1 — e, for the 
global-addressing technique presented here. This higher 
fidelity requires larger pulse duration for a fixed peak 
Rabi frequency. 

The characteristic pulse width for a repeated-local- 
addressing technique with a single ion requires A'^ steps 
with an N times lower error [cf. Eq. (PT]) ]. 

^ ^ V2NHN/e) ^ VN HN/e) 
TTV^Q y2 7r.g(0) 

The total interaction time for a sech pulse is again 
Tlocai ~ lOT. The single-ion implementation therefore re- 
quires longer time than the global-addressing technique, 
by a factor of ln(A^/e)/ln(l/e). For example, for e = 0.01 
and = 8 ions, this factor is about 1.5. Therefore the 
number of heating events during the preparation of the 
Fock states increases by the same factor. 

In summary, our global-addressing technique allows for 
a shorter time duration than single-ion techniques, and 
in addition, offers an appealing simplicity of implementa- 
tion, with the need to engineer a single interaction only. 

Finally, we consider the preparation of motional states 
of a single ion, as shown in Fig. [7Ji, whereby a sequence of 
adiabatic pulses can be applied to add or remove phonons 
one at a time. In this respect the proposal is similar to 
earlier work involving single ions, such as [l6j which used 
sequential pulses of precise area. However, our method 
possesses a significant advantage over previous work - 
namely a natural robustness with respect to experimen- 
tal imperfections. This robustness arises from the adia- 
baticity of the process, as discussed in Sec. lIIIl and should 
allow a significant increase in the fidelity of the prepara- 



tion procedure. 



VI. CONCLUSIONS 

We have shown how vibrational Fock states of a 
trapped-ion chain can be produced in a single step and 
using a single laser pulse. This is achieved by sweeping 
the laser frequency through resonance with a motional- 
sideband of the centre-of-mass vibrational mode. The 
key theoretical ingredient of our proposal is a particu- 
lar choice of basis states in which to analyse the problem 
(the generalized Morris-Shore basis), wherein the Hilbert 
space factorizes into a collection of independent ladders. 
Sweeping the laser detuning through resonance induces 
a 'bow-tie' multiple-level crossing. With the proviso that 
the pulse and detuning satisfy an adiabaticity condi- 
tion, this pulse connects the lowest and highest levels of 
each ladder. Moreover, the zero-phonon and A^-phonon 
Fock states arc located at the opposite ends of the longest 
chain, which provides an extremely simple recipe for cre- 
ating phonon Fock states in a single step. Since the state 
preparation takes place adiabatically, it is naturally ro- 
bust against fluctuations in the experimental parameters, 
such as the laser detuning and chirp, the Rabi frequency, 
and the pulse duration. Numerical simulations show that 
the technique also performs well on the boundary of the 
Lamb-Dickc regime. 

Our technique may also be used to prepare GHZ states 
of the ions' internal degrees of freedom and superposi- 
tions of motional Fock states of the whole chain. We 
estimate that the fidelity of each step could be over 99% 
for an eight-phonon Fock state if the centre-of-mass mode 
and a single laser pulse are used, and we have also sug- 
gested several ways in which heating effects could be fur- 
ther reduced. 
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These are the 0{rf') terms; there are also higher-order 
terms, which are smaller again in magnitude. 
The rotating-wave approximation is easy to satisfy in 
the adiabatic limit, since any transitions driven off- 
resonantly on the carrier transition cancel out over the 
course of the pulse. Therefore, Q, < u/\/N is enough to 
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